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Identification  of  Free  Oscillation  Spectral 
Peaks  Using  the  Excitation  Criterion 

In  this  final  annual  report,  we  describe  the  work 
of  Jorge  Mendiguren  on  the  identification  of  free  oscillation 
spectral  peaks  by  the  new  powerful  method  called  "excitation 
criterion".  An  intensive  search  of  higher  modes  was  done 
using  data  for  the  July  30,  1970  Colombian  deep  shock. 
Spheroidal  oscillations ,  ranging  from  the  fundamental  to 
the  sixth  higher  overtone,  and  torsional  oscillations  up  to 
the  third  higher  overtone  were  identified.  Spheroidal 
oscillations  for  periods  smaller  than  100  seconds  were  clearly 
identified  for  the  third  and  fifth  overtone  . 

The  standard  deviation  of  the  measurements  was  estimated 
applying  the  excitation  criterion  to  different  subsets  of 
stations .  Some  of  the  eigenperiods  measured  for  higher 
overtones,  especially  the  third  and  fifth  overtones  at 
periods  near  100  seconds  were  found  to  be  very  stable.  Dis¬ 
crepancies  between  eigenperiods  for  different  subsets  of 
stations  are  as  small  as  0.01%.  This  is  in  agreement  with 
the  idea  of  a  lower  mantle  with  less  important  lateral 
heterogeneities  than  the  upper  mantle. 

Comparison  of  the  eigenperiods  determined  for  both  the 
Maskan  1964  and  the  Colombian  1970  earthquakes  showed 
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systematic  differences.  That  difference  is  of  the  order  of 
0.3*  for  oT  modes  while  for  qS  modes  it  is  much  smaller, 
around  0.04*.  There  are  indications  that  such  a  discrepancy 
IS  also  present  for  overtones.  This  means  that  the  lateral 
heterogeneities  of  the  earth  are  such  that,  even  with  good 
coverage  of  stations,  the  average  eigenperiods  measured  for 
shocks  with  different  epicenters  and  source  mechanisms  will 
be  different  by  a  measurable  amount. 

The  effect  of  the  source  finiteness  on  the  spectral 
resolution  was  theoretically  analyzed.  We  found  that  for 
oscillations  with  wavelengths  about  6  times  the  source  linear 
dimension  or  longer,  the  effect  of  source  dimension  is  not 
too  important.  A  simple  point  source  can  be  assumed  if  we 
discard  those  stations  vhere  the  theoretical  amplitude  of 
the  mode  under  study  is  ~2.5  times  smaller  than  the  average 
amplitude  computed  for  all  the  stations. 

The  follov/ing  paper  by  Mendiguren  is  reprinted  from  the 

Sg2B*gsical  Journal  of  the  Royal  Astronomical  Society,  33, 
1973. 
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Ctophys.  1.  R-  astr.  Soc.  (1973)  33,  281-321. 


Identification  of  Free  Oscillation  Spectra.  Peaks  for 
1970  July  31,  Colombian  Deep  Shock  usins  the 
Excitation  Criterion 

Jorge  A.  Mcndiguren 
(Received  1973  February  26)* 


Summary 

The  excitation  criterion  was  applied  to  identify  spectra1  peaks  of  free 
oscillations  excited  by  the  Colombian  deep  .took  ol  1  )70  July  a I. ^ 

160  new  overtones  were  identified.  They  ranee  from  ,  f  to  3T  tors  o na 
modes  and  from  ,S  to  hS  spheroidal  modes.  Some  ot  the  mooes  iden  died 
have  periods  below  90s.  The  standard  deviation  ol  the  measured  e  _.ei  - 
ncriods  is  in  eeneral  smaller  titan  0-1%.  Comparison  ot  ,,S  and  0T 
eigen  periods  measured  for  the  Alaskan  (1964)  and  Colombian  shocks 
shows  a  systematic  difference.  The  difference  is  larger  than  the  error  ot  the 
measurement,  and  this  indicates  that  the  average  ementmquenc.es  measured 
for  each  event  are  biased  by  the  epicentre  location.  The  elicit  ot  lateral 
heterogeneities  of  the  Earth,  aftershocks,  and  the  choice  ol  source  and 
Earth  model  when  using  the  excitation  criterion  are  discussed  in  detail 


1.  Introduction 

In  recent  year,  there  ha,  been  an  increasin'  emphasis  on  the  u« icf  free 
data  to  derive  models  of  the  Earth's  interior,  e  g.  Press  (1968),  Haddon  and  bum 
(1969),  Dziewonski  &  Gilbert  (1972a).  Johnson  ( 1972).  Most  ol  thc  ‘(J* L^O-Sos 

data  available  for ’hose  studies  correspond  to  modes  with  periods  longer  tha  -  -  ■ 

BM  itis  clear  that  additional  h, alter  mode  data  wtll  be  necessary  tor  a  more  detaded 
determination  olthe  Earth  strucldrc  tWteeins  W2>.  and  theretore  net*  tdentihca 

nf  overtones  at  shorter  periods  are  highly  desirable. 

Unfonunatcly  for  periods  below  250-300  s  it  is  not  poss.ble  to  obtain  unambitious 

£fe,  r =1 

Slichter  1961),  or  separation  ol  shear  and  volumetric  strains  te._.  t 
c—uu  1966).  Some  of  the  limitations  of  these  methods  can  be  avoided  using 
differential  attenuation  criterion  proposed  bv  Dratler  et  al.  ( '971).  an  ^n^lve  y 
used  by  Dziewonski  &  Gilbert  ( 1972b)  lor  periods  below  250  s.  This  method  s  . 
on  the  difference  in  attenuation  for  different  modes,  and  can  be  applied  idinti  . 

at  Short  periods  is  the  criterion  based  on  tlte 
ray-mod;  duality;  originally  proposed  by  Brune  (1964),  it  was  later  used  by  Als  p 

•Rece'  /ed  in  original  form  1973  February  9. 
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Brune  (1965)  and  Nowroozi  (1972).  The  identifications  obtained  using  this  method 
have  not  been  critically  tested  yet,  and  additional  applications  are  necessary  for  a  clear 
evaluation  of  its  practical  merits. 

Unlike  the  earlier  methods,  the  excitation  criterion  (Mendiguren  1972a,  b,  c,  1973) 
renders  reliable  identifications  at  all  periods  and  can  be  applied  to  torsional  as  well  as 
spheroidal  modes,  even  in  those  cases  with  low  Q  values.  Theoretically  this  method 
should  work  for  any  well-excited  mode  if  its  amplitude  and  phase  at  the  observing 
stations  can  be  predicted  using  a  suitable  Earth  and  source  model. 

In  the  following  sections,  the  excitation  criterion  is  discussed  in  more  detail  than 
in  former  publications.  And  finally  the  results  of  its  application  to  the  free  o  cillations 
excited  by  the  Colombian  deep  shock  of  1970  July  31  are  thoroughly  displayed. 


2.  Theoretical  and  observed  spectra  of  free  oscillations 

The  basic  assumption  underlying  the  excitation  criterion  is  that  the  observed 
spectra  of  free  oscillations  can  be  predicted  using  suitable  models  for  the  Earth  and  the 
source  mechanism. 

In  this  section  it  will  be  shown  that  for  the  case  of  the  free  oscillations  excited  by 
the  Colombian  deep  shock  of  1970  July  31.  a  simple  dcuble-couple  source  in  a  laterally 
homogeneous  earth  is  an  adequate  model  to  predict  the  spectra  of  the  fundamental 
spheroidal  modes.  In  later  sections  this  simple  model  is  successfully  used  in  the 
excitation  criterion  scheme  to  identify  higher  modes  uniquely,  and  this  indirectly 
proves  that  this  model  is  also  suitable  to  predict  the  spectra  of  overtones.  Saito's 
formulae  (Saito  19o7)  have  been  used  in  this  study  to  compute  the  theoretical  spectra. 

The  displacement  spectra  of  free  oscillations  excited  by  a  double-couple  source, 
having  a  step  function  as  source  time  function  and  being  imbedded  in  a  laterally 
homogeneous  earth  can  be  expressed  as  follows  (Saito  1967): 


(ir  =  I  I  I— 7,T(a)  T”  (0,<l>)  ela"-  S 

«  m-0  J  a>  „  J 


[is  =  IZ  Z— 

m  mm0  J  0}  „  j 


{  7,s(a)  S,"  .  (0,  <t>)  +  7/  (a)  S2m,  M  <f>)} 


(1) 

(2) 


where 


r,0,4>  «  spherical  co-ordinate  system:  r  is  measured  from  the  centre  of  the 
Earth;  the  0  =  0  axes  passes  through  the  hypoccntrc;  <p  is  measured 
counterclockwise  from  the  fault  strike  direction. 

j  -  radial  order  number. 

n  =  co-latitudinal  order  number;  order  of  Legendre  polynomial. 
m  =  azimuthal  order  number;  degree  of  Legendre  polynomial. 
a  =  Earth  radius. 

( o,j  «  angular  eigenfrequency  corresponding  to  the  mode  of  radial  order  j. 

t  =  time  measured  from  the  source  origin  time. 

S,T  =  superscripts  indicating  spheroidal  and  torsional  modes. 

7,*,  Tj1  =  radial  and  0  displacement  eigenfunctions  for  spheroidal  modes. 

7,T  =  (p  displacement  eigenfunction  for  torsional  modes. 

A,  B  =  coefficients  which  are  functions  of  focal  depth,  Earth  model  and 
orientation  of  the  double  couple. 

T,m  (9,  tf>),  S  ",  „  (0,  <p),  Si"  „  (0, 4>)  =  vector  spherical  harmonics 


Identification  of  fret  oscillation  spectral  peaks 
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4> 

1  3  7- 

-ay- 

sinfr  dip 

36 

0 

0 

ay- 

1  ay„- 

36 

sin 0 .  dtp 
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/COS  , 

Y„m(0,<i>)=  P„"(cos0)  (  m4>) 

'sin  ' 


The  hypocentre  parameters  used  to  compute  1  and  2  are  the  se  given  by  NOAA 
Origin  time  17h08m5-4s 

Latitude  1-5°  S 

Longitude  72-6°  W 

Depth  651  km 

(Magnitude  n,b  =  7-1) 


N 


Flo.  I.  Equal  area  projection  of  the  lower  hemisphere.  Open  and  full  circles 
indicate  dilatations  and  compressions,  respectively. 
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«k.Hie.r,0dal/fJan!  S0lu,i0n  nceded  10  comPu,e  A  and  B  was  derived  from  P  wave 
observations  (Mendigurcn  1972a)  and  it  is  shown  in  Fig.  1.  Its  parameters  are: 


Dip  direction 

Dip 

Trend 

Plunge 

Plane  a 

N 122° W 

58° 

Plane  b 

N75°E 

32° 

Pressure  axis 

N32°E 

75° 

Tension  axis 

N116°W 

12° 

Null  axis 

N153°E 

8° 

J.h' HBI  '6.Earcth  model  of  Haddon  &  Bullen  (1969)  was  adopted  in  this  study  but 

6-  the  Select!on  the  Earth  model  is  not  critical  in  the  excitation 
criterion  scheme.  The  corresponding  free  oscillations  eigenfunctions  were  computed 

for  M  W rtUn  V?  Wi?si,ii  Theoretical  spectra  were  commuted 

s'atl0ns,  wh,le  ,he  observed  snec  ra  was  obtained  from  long  period 

S  31  th°Se  s,a,l0"s-  A  l,st  of  the  stations  is  given  in  the  Appendix. 
I  he  method  of  data  processing  is  also  described  in  the  Appendix. 

2 if"d  •jS,h0w  'he  rcsuUinP  theoretical  radial  displacement  spectra  forfunda- 
cntal  spheroidal  modes  compared  to  the  observed  ones  at  20  stations.  The  main 
features  of  the  observed  spectra  are  closely  reproduced  by  the  theoretical  spectra  of 
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^iiJilLiLIJlUlLl.liliil,  Lid--.,, 

Fw.  2  Continuous  Hnes  indicate  the  observed  radial  displacement  spectrum, 
ertical  bars  arc  tne  theoretical  spectral  lines.  A  is  the  epiccntra  I  distance.  Azisthe 
azimuth  at  the  epicentre  measured  clockwise  from  nonh.  The  patter  of  alternating 
large  anu  small  amplitude  peaks  is  characteristic  of  stations  near  A  *  90°. 
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Mratiflcstioo  of  free  oscillation  spectral  peaks 
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Al' 


Flo.  3.  Continuous  lines  indicate  the  observed  radial  displacement  spectrum. 
Vertical  bars  are  the  theoretical  spectral  lines.  is  the  epicentral  distance.  Az  is 
the  azimuth  of  theepicentrc  measured  clockw  isc  from  north.  Note  the  shift  in  nodes 
between  the  (I  and  R  spectra  for  KOD. 
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STATION  LAH 
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Fig.  4.  Theoretical  line  spectra  of  the  radial  displacement  of  spheroidal  oscillations 
at  LAH.  Numbers  under  each  line  indicate  the  co-laiiiudinal  order  number. 
+  and  1  indicate  that  ihc  phase  of  the  mode  is  0  or  tt  respectively. 


oS  modes.  The  QS  spectra  explain  by  itself  most  of  the  observed  spectral  features 
because  as  shown  in  Fig.  4  for  LAH  station,  nS  dominates  the  vertical  displacement 
at  that  frequency  range.  The  amplitude  of  the  other  modes  is  much  smaller.  The  same 
occurs  at  all  other  stations 

The  close  agreement  shown  in  Figs  2  and  3  indicates  that  a  point  source  in  a 
laterally  homogeneous  Earth  is  an  adequate  model  to  predict  the  actual  displacement 
spectra  for  this  e\ cnt.  Limitations  of  this  simple  model  due  to  lateral  heterogeneities 
of  the  real  Earth  and  source  finiteness  will  be  discussed  in  later  sections. 

3.  The  excitation  criterion  applied  to  the  spectra  observed  at  a  single  station 

The  close  agreement  shown  in  Figs  2  and  3  suggests  that  a  simple  method  to 
identify  spectral  peaks  is  to  compare  theoretical  to  observed  spectra  at  single  stations. 
This  technique  was  applied  to  »he  oscillations  excited  by  the  Colombian  shock  and 
many  overtones  were  identified.  Some  examples  of  those  identifications  on  spectra 
obtained  from  records  of  NVWSSN  stations  and  the  NE  strainmeter  at  Nana  (Peru) 
are  shown  in  this  section. 
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Fra.  5.  Theoretical  and  observed  spectra  of  the  radial  displacement  component  at  KEV 
The  numbers  under  the  theoretical  line  spectra  indicate  the  co-latitudiml  order  number. 
*+*  and  *  I  ’  indicate  that  the  phase  of  the  mode  is  0  or  n  respectively. 


Fig.  5  shows  part  of  the  radial  component  spectrum  observed  at  KEV  compared 
to  the  theoretical  spectra  o‘- different  spheroidal  modes.  The  identification)  of  the  ,S 
peaks  are  based  on  the  dominant  excitation  expected  for  those  modes  ard  the  agree¬ 
ment  between  the  theoretical  shape  of  the  ,S  spectrum  and  the  main  features  of  the 
observed  one.  The  other  modes  arc  not  so  well  excited  and  the  spacing  between 
theoretical  peaks  does  not  match  the  observed  one. 

In  a  similar  way.  Fig.  6  shows  ;S  peaks  identified  in  the  co-latitudinal  displacement 
spectra  at  ESK.  5S  modes  identified  at  ATU  are  shown  tn  Fig.  7. 

Usually  the  ,S  peaks  in  the  colatitudtnal  component  spectra  are  confusingly  mixed 
with  0S  peaks  in  the  frequency  range  shown  tn  Fig.  8,  but  in  the  particular  case  of 
BOG  the  extended  node  in  the  0S  spectra  makes  the  3S  peaks  clearly  visible.  The 
theoretical  prediction  of  the  node  in  „S  spectra  gives  support  to  those  identifications. 
At  both  sides  of  the  node  the  identifications  are  iess  reliable. 
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Fig.  6.  Theoretical  and  observed  specira  of  the  co-laiitudinal  displacement 

component  at  ESK. 


The  spheroidal  modes  have  an  azimuthal  displacement  component  of  very  small 
amplitude  when  compared  to  the  radial  orco-latiludinal  ones  and  their  corresponding 
peaks  nave  never  been  observed  before.  As  their  amplitude  goes  like  I  'sin  f?  (see  the 
<jy  component  of  S"\  „  in  expression  |3)  they  only  eould  be  observed  close  to  the 
epicentre  or  its  antipodes.  Fig.  y  shows  some  of  those  peaks  eorresponding  to  the 
0S  modes  identified  in  an  extended  node  of  the  0T  spectra  at  BOG,  onK  6t.1  km  from 
the  epicentre.  At  both  sides  ol  the  node  of  „  T  the  peaks  are  confusedly  intermixed 
and  no  clear  identification  is  possible.  Fig.  10  shows  a  more  complete  theoretical 
spectra  of  the  azimuthal  component  at  BOG.  Only  the  knowledge  of  the  theoretical 
spectra  makes  those  identifications  reliable. 

The  strain  spectra  observed  by  the  NE  component  strainmeter  at  NNA  is  shown 
in  Fig.  1 1  in  comparison  with  the  theoretical  one.  As  the  NE  component  strainmeter 
is  directed  at  an  angle  a  -  17  from  the  station  to  epicentre  direction,  its  spectrum 
shows  intermixed  torsional  and  spheroidal  peaks.  The  theoretical  strain  was  computed 
as  a  function  of  the  cou ,  c0+,  c ++  strains  as  follows. 

e„  =  ctl  cos2  a  +  C++  sin2  a  4-  ce+  sin  a  cos  a 

eH ,  et<l  and  C++  were  computed  taking  appropriate  derivatives  of  the  displacement 
given  by  expressions  (I )  and  (2).  The  agreement  shown  in  Fig.  1 1  indicates  that  our 
simplified  source  and  Earth  models  arc  also  valid  at  the  lowest  frequency  range. 
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Fra.  8.  Theoretical  and  observed  spectra  of  the  co-latitudinal  displacemen 

component  at  BOG. 
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Fig.  10.  Theoretical  spectra  of  the  azimuthal  displacement  component  at  BOG. 
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Fto.  U  Observe!  and  theoretical  strain  spectrum  lor  the  NMV  E  strammeter  at 

T,  r°bSerI  5pean,m  ,s  lhc  secoml  fr™  top  The  spectrum  at  .he 
top  results  from  the  convolution  ot  the  theoretical  line  spectrum  v.uh  a  sin  r  v 

.r2?.l|°naSlnlUlaT8  'he  e,,cc'  of  ,hc  rtni,c  record  length  nT  modes  luve  been 
slightly  displaced  touard  lower  frequencies  to  match  the  observations  „Sn  and 
»■>,  to  require  larger  amplitudes  to  match  the  observations.  Those  observed 
peaks  indicated  with  a  question  mark  have  not  a  theoretical  equivalent  and  are 
attributed  to  recording  noise. 
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All  these  examples  indicate  that  if  the  spectra  at  only  a  few  stations  are  available, 
some  of  the  observed  peaks  can  be  identified  by  comparing  theoretical  to  observed 
amplitude  spectra  at  sinclc  stations.  Lisuallv  this  approach  renders  more  convincing 
identifications  when  only  one  mode  dominates  the  spectrum,  like  the  2S  modes  in 
Fig.  6.  Those  modes  with  smaller  amplitude  cannot  be  unambiguously  identified 
using  this  technique.  In  those  cases  the  method  described  in  the  next  section  should 
be  used. 


4.  The  excitation  criterion  at  plied  to  spectra  observed  at  many  stations 

The  spectra  indicated  as  '.ir|,  M„|  and  |/1*|  in  Figs  12  and  13  are  the  absolute 
value  of  radial,  co-latitudinal  and  azimuthal  spectra  respectively,  summed  over  all 
stations.  This  summation  does  not  help  much  to  identify  particular  spectral  peaks 
but  enhances  modes  which  arc  well  excited,  like  those  peaks  between  0  02  and 
0-055  cycles/sec  in  |,4r|  or  the  sequence  of  6S  peaks  in  the  same  spectrum  between 
0-07  and  0-008  cycles  see.  Similarly  a  clear  sequence  of  equally  spaced  5S  peaks  is 
observed  in  |/lf)|  between  0-0075  and  0-0086 cycles/scc. 

A,  in  Figs  12  and  13  is  the  spectrum  resulting  from  vectorial  summation  (phase 
and  amplitude  taken  into  account)  of  the  same  radial  spectra  used  in  |.4,|.  The 
difference  between  A,  and  |^rl  is  striking,  most  of  the  peaks  well  developed  in \Ar\  have 
disappeared  in  Ar.  This  is  due  to  the  fact  that  for  t  =  0  the  phase  of  a  peak  at  a  given 
station  can  be  0  or  r  (see  expression  (1)  and  (2)).  The  phase  of  a  given  peak  is  n  at 
some  stations  and  0  at  others  and  when  the  spectra  are  vectoriallv  added  those  peaks 
tend  to  cancel  each  other.  0S:,  in  Ar  is  a  clear  exception;  it  shows  a  large  amplitude 
because  it  just  happens  that  it  has  the  same  phase  at  almost  all  the  stations. 

It  is  possible  to  make  the  phase  of  a  given  mode  the  same  in  all  the  spectra  by 
adding  j:  to  the  observed  jpcctra  of  those  stations  where  its  theoretical  phase  is  n  and 
leaving  the  remaining  spectra  unchanged.  By  doing  that,  the  phase  ot  the  mode 
becomes  equal  to  0  at  all  stations  and  its  peaks  will  add  constructively  when  vcctorially 
added.  The  result  of  such  a  summation  after  phase  shift  lor  the  0S:o  mode  is  shown 
as  Ar*P  in  Figs  12  and  13.  Now  0S:0  dominates  the  spectrum;  most  of  the  other 
peaks  have  almost  disappeared  because  they  have  a  difierent  phase  at  different 
stations.  Note  that  0S:o  has  disappeared  now.  Some  peaks,  such  as  0S:2  at  the  right 
of  0S20,  have  intermediate  amplitudes  because  they  have  a  spatial  phase  distribution 
close  to  that  of  C,S:0.  Then  it  is  clear  that  having  a  large  number  of  spectra  of  stations 
all  over  the  world  and  finding  a  suitable  source  and  Earth  model  to  predict  the  spectra 
at  each  station  it  is  possible  to  filter  out  and  unambiguously  identify  any  peak  if  it  is 
sufficiently  well  excited.  Figs  14  to  28  show  the  result  of  application  of  this  method 
to  cases  where  the  author  believes  that  the  identification  is  correct  and  the  measured 
eigenfrequenev  is  accurate.  Doubtful  cases  arc  not  shown.  Particular  comments  on 
each  mode  identification  are  given  in  the  Appendix.  The  display  ot  detailed  results 
has  the  purpose  of  allowing  a  critical  analysis  of  these  identifications  when  used  for 
inversion. 

Table  1  gives  the  eigenperiod  measured  in  each  case.  The  cigenfrequency  was 
found  by  polynomial  interpolation  over  4  points  of  the  discrete  amplitude  spectra. 
In  order  to  enhance  the  signal  to  noise  ratio  the  spectra  at  those  stations  where  the 
theoretical  amplitude  of  the  mode  under  study  was  less  than  0-5  times  the  average 
amplitude  were  discarded  before  summation.  The  resulting  eigcnperiods  arc  an 
amplitude  weighted  average  over  all  observing  stations.  This  scheme  has  the  effect  of 
assigning  weights  proportional  to  the  signal  to  noise  ratio.  Identifications  only  up  to 
the  sixth  spheroidal  overtone  and  third  torsional  overtone  were  attempted  in  this  study. 


frequency  ( cycles/secono ) 

Flo.  12.  Spectra  resulting  from  summation  of  spectra  observed  at  many  stations 
over  the  Earth.  The  continuation  of  these  spectra  for  smaller  periods  is  shown  in 

Fig.  13. 
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Table  1 


Eigenperiods  found  using  the  excitation  criterion.  *  T'  and  '  S  D.'  are  the  measured 
eigenperious  and  the  corresponding  standard  deviations  in  seconds.  '  Delta  '  is  the 
difference  in  seconds  between  the  eigenpertod  measured  using  all  stations  and  the  average 

over  four  subsets  of  stations. 


Observed  eigenperiods 


oS. 


n 

T 

S.D. 

Delia 

13 

472 • 979 

0-072 

-0-000 

14 

448 '201 

0  074 

0-013 

13 

426  056 

0  050 

0-000 

16 

406' 746 

0-103 

-0-022 

17 

389  134 

0-118 

0  007 

18 

374  010 

0  048 

0-042 

19 

360112 

0-056 

—  0-011 

20 

347 '505 

0  057 

0-010 

21 

335-811 

0-100 

-0  007 

22 

325  059 

0  027 

-0  004 

23 

315-213 

0  116 

-0  021 

24 

306-279 

0-065 

0  030 

23 

297-664 

0-013 

0  005 

26 

289-601 

0  026 

-0  000 

27 

282-181 

0  054 

0  005 

28 

275-113 

0  054 

-0  024 

29 

268-436 

0  043 

0  009 

30 

262-060 

0  098 

0  013 

31 

255-952 

0-110 

0-009 

32 

250-309 

0-095 

0  019 

33 

244-919 

0-133 

0-030 

34 

239-585 

0  041 

0  012 

33 

234-579 

0  036 

0  008 

36 

229-806 

0  038 

-0  008 

37 

225-215 

0-036 

-0  009 

38 

220-742 

0015 

-0-002 

39 

216-478 

0  033 

0-001 

40 

212-381 

0-035 

-0  004 

41 

208-347 

0  029 

0-003 

42 

204-579 

0  022 

0-001 

43 

200-927 

0-031 

0-004 

44 

197-398 

0-028 

-0-000 

43 

193-875 

0-041 

-0  002 

46 

190-624 

0  042 

0-015 

47 

187-258 

0  093 

-0-040 

48 

184-289 

0139 

0-052 

49 

181  002 

0-021 

0  004 

30 

178-351 

0-115 

0-037 

31 

175-267 

0-075 

-0-017 

32 

172-541 

0-032 

0  003 

33 

170-089 

0-011 

-0-002 

34 

167-368 

0-047 

0-005 

33 

164-840 

0-064 

0-015 

36 

162-482 

0-130 

0-033 

57 

160-265 

0-066 

0-076 

38 

157-844 

0-060 

-0-006 

39 

155-011 

0  032 

-0-002 

>5. 


II 

T 

S.D. 

Delta 

16 

298-930 

0-398 

-0-047 

18 

274-752 

0-282 

0  035 

21 

245-017 

0-278 

0-004 

.5. 


n 

T  . 

S.D. 

Delta 

23 

229-347 

0-057 

0  006 

29 

190-894 

0-115 

0  000 

36 

161-353 

0  087 

-0-058 

38 

154-758 

0-079 

0  040 

39 

151-475 

0-001 

-0-001 

41 

145  832 

0  052 

0-007 

42 

143-047 

0-121 

—  0  103 

44 

138-559 

0  410 

-0-048 

50 

125  386 

0-282 

—  0  118 

52 

121-955 

0  032 

-0  013 

53 

120-074 

0-034 

0  014 

55 

116-581 

0  052 

0  009 

iS, 


n 

r 

S.D. 

Della 

6 

595-011 

0-756 

0  039 

7 

535-260 

0  990 

0-279 

10 

415-916 

0  770 

-0-279 

II 

387-999 

0  207 

-0-082 

12 

365-230 

0-248 

0-009 

13 

343-219 

0-357 

-0-170 

14 

325-867 

0-462 

-0-042 

15 

306  400 

0-407 

-0  336 

28 

169-251 

0  033 

-0  002 

29 

164-587 

0-097 

-0-043 

30 

160-553 

0  031 

0-010 

31 

156-616 

0-094 

-0-043 

32 

152-684 

0-037 

0-010 

33 

149-169 

0-155 

0-094 

35 

142-614 

0-074 

-0-008 

36 

139-507 

0-041 

0-002 

37 

136-648 

0-001 

-0  000 

38 

133-858 

0-104 

0  086 

39 

131-134 

0-060 

0-016 

40 

128-538 

0-002 

-0  001 

sS. 


n 

T 

S.D. 

Delta 

7 

371-856 

0  108 

-0  075 

8 

354-296 

0183 

0  021 

9 

338-464 

0-803 

0-374 

10 

323-954 

0-152 

-0-057 

11 

310-149 

0-260 

0-174 

12 

297-475 

0-075 

-0  006 

13 

284-982 

0-050 

-0  029 

14 

273-351 

0-098 

-0  007 

15 

262-467 

0  091 

0  031 

16 

251-984 

0  043 

-0  003 

17 

242-443 

0-125 

-0-024 

18 

233-286 

0  120 

-0-077 

19 

224-910 

0-121 

-0-019 

-20- 
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Table  1  ( continued ) 


H 

T 

»s. 

S.D 

Delta 

n 

T 

S.D. 

Delta 

20 

21 

22 

23 

24 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

216-949 
209-564 
202-620 
195  938 
190  068 
111-451 
109  476 
107-776 
105-976 
104  197 
102-585 
101  141 
99-462 
98-014 
96-439 
94-961 

0-064 

0-035 

0-047 

0  295 

0  033 
0-033 
0-082 
0-032 
0-094 

0  101 
0-074 

0  032 

0  084 
0-020 
0-069 
0-052 

0-010 

0  014 

0  001 
0-066 
-0  011 
-0-001 

0  039 

0  010 
0-049 
-0-070 
0-029 
0-008 
0-059 
-0-006 
-0  013 
-0-036 

16 

17 

18 

19 

in 

21 

22 

23 

24 

25 

26 

27 

28 

29 

32 

172-342 
166-159 
160-577 
154-973 
! 50  090 
145-756 
141-909 
138-329 
134-794 
131-784 
128-697 
125-951 
123-514 
121-043 
114-452 

0-104 

0  101 
0-218 
0-084 
0-045 
0-017 
0-031 
0-053 
0-053 
0-034 

0  101 
0-036 
0-080 
0-019 
0-073 

-0-079 
-0  026 
-0-183 
-0-063 
0  030 
0-000 
-0-007 
0  049 
-0-020 
-0-018 
-0-072 
0  007 
0  001 
-0-015 
0-032 

53 

93-725 

0-193 

0-161 

54 

92-336 

0  184 

-0-013 

oT, 

56 

90-103 

0  001 

-0-000 

n 

T 

S.D. 

Delta 

57 

88-882 

0-000 

0-000 

11 

574-487 

0  181 

0-041 

58 

87-654 

0-031 

-0  011 

12 

539-946 

1  233 

0-155 

13 

506-832 

0-306 

0-183 

.  c 

14 

477-511 

0-281 

0  141 

n 

T 

S.D. 

Delta 

15 

16 

452-933 

431-229 

0-220 

0-590 

-0  018 
-0-070 

32 

135-649 

0-071 

-0-001 

17 

410-350 

0  210 

-0  095 

33 

132-825 

0-102 

0-011 

18 

392-354 

0-259 

0  098 

34 

129-931 

0  063 

-0  014 

19 

374-555 

1-086 

-0-738 

35 

127-243 

0-030 

0  007 

20 

358  718 

0  174 

-0-033 

36 

124-668 

0  045 

0-006 

21 

346- 142 

0-286 

-0  028 

37 

122-242 

0  019 

-0-001 

22 

332-855 

0-251 

0  110 

38 

119-843 

0-035 

0-009 

23 

321-953 

0  351 

0-002 

39 

117-608 

0-035 

0-025 

24 

310-362 

0-322 

-0-077 

40 

115-436 

0-077 

0  019 

25 

300186 

0-260 

-0-022 

41 

113-375 

0-021 

0  001 

26 

290-660 

0  079 

-0-047 

27 

281-887 

0-321 

0-198 

29 

264-934 

0-178 

-0-024 

30 

257-830 

0-325 

0  083 

n 

T 

S.D. 

Delta 

5 

e 

E 

s 


s 

f. 

2 

s 

s 

l 


24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 
46 


147  027 
143-594 
140-266 
137-122 
134  064 
131-176 
128-477 
125-881 
123-433 
121-030 
118-768 
116-633 
114-527 
112 543 
110-590 
108  665 
106-850 
105-051 
103-354 
101-707 
100-079 
97-040 


0-028 
0  019 
0-066 
0-047 
0-057 
0-020 
0  061 
0-056 
0-053 
0-017 
0-020 
0-040 
0-038 
0-017 
0  000 
0-001 
0-001 
0  026 
0  002 
0-052 
0-001 
0-016 


0  005 
0  004 
0-008 
0-023 
0  005 
0  010 
0  030 
0  011 
026 
006 
000 
006 
002 
009 
000 
000 
000 
015 
000 
032 
000 
009 


.r. 


S.D. 


Dtlta 


7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 


475-191 
438-962 
407-752 
381-228 
359-077 
339  239 
322-837 
306-996 
293  364 
280-515 
269-202 
259-092 
249-158 
240-977 
232-529 
225-190 
218-379 
211-948 
205-853 
200-270 


0- 

0- 

0 

0 

0- 

0 

0- 

0- 

0- 

0- 

0- 

0- 

0- 

0- 

0- 

0 


•617 
•393 
•797 
•504 
•486 
151 
•398 
•197 
•074 
151 
•253 
•343 
155 
•224 
•112 
052 
0  141 
0-024 
0-103 
0-040 


-0- 

0- 

-0- 

0- 

-0 

0- 

-0 

0- 

0 

-0 

0- 

0 

0- 

0 

-0- 


■096 
•067 
■307 
■098 
•006 
•172 
•075 
•009 
CC6 
•054 
113 
•138 
■044 
•048 
•001 
-0  010 
-0  011 
0  007 
0-034 
0-001 
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Table  1  ( continued) 

«r.  _  ,r. 


It 

T 

S.D. 

Delta 

n 

T 

S.D. 

27 

194-943 

0-022 

-0-016 

25 

17* • 123 

0-126 

28 

190  085 

0-109 

-0-035 

27 

162-581 

0-138 

29 

185-338 

0-096 

-0-020 

28 

158-409 

0-150 

31 

176-849 

0-130 

0  053 

29 

154-571 

0-049 

32 

173-027 

0  324 

-0-205 

32 

144-201 

0-210 

33 

169  156 

0-041 

0  019 

37 

130-406 

0-  268 

34 

165-715 

0-085 

0-064 

38 

128-168 

0-099 

35 

162-356 

0-032 

-0-008 

41 

121-571 

0-048 

36 

159  113 

0  022 

0-018 

42 

119-330 

0-172 

37 

155-982 

0  034 

0-013 

38 

153-165 

0-115 

-0-060 

39 

150-253 

0  056 

0-00S 

40 

147-681 

0-007 

0-007 

%T. 

SJ>. 

41 

145-121 

0  049 

0-006 

H 

T 

42 

142-661 

0-052 

0  010 

ft 

4 

43 

140-230 

0-111 

0  006 

17 

189-973 

0-247 

44 

137-960 

0-0S2 

0-021 

19 

178-104 

0-080 

45 

135-638 

0-330 

-0-285 

24 

154-809 

0-184 

47 

131-592 

0-222 

-0-101 

25 

150-663 

0-047 

49 

127-651 

0  021 

-0  012 

26 

146-936 

0-310 

51 

124-127 

0-539 

-0-411 

28 

140-390 

0  099 

54 

119  127 

0-050 

-0-004 

29 

137-193 

0-067 

299 


Delta 

-0  033 
0  000 
-0  037 
0-005 
0016 
-0172 
0  055 
-0  007 
0-177 


Delta 

-0-154 
-0  016 
-0-110 
-0-025 
-0  189 
-0-067 
0-020 


5.  Accuracy  of  the  measured  eigenperiods 

The  error  in  the  eigenperiods  determined  by  vectorial  summation  of  spectra  was 
estimated  as  follows.  The  stations  were  grouped  into  four  subsets  based  on  the 
geographical  location  with  respect  to  the  epicentre.  The  lirst  subset  includes  those 
stations  located  at  epicentral  distance  between  0  and  5”.  10'  and  15’.  20a  and  25'  and 
so  on.  The  second  subset  includes  the  remaining  stations.  The  third  subset  was 
formed  with  those  stations  located  at  azimuth  between  0  and  5°.  10  and  15’,  20’  and 
25°  and  so  on.  The  fourth  subset  includes  all  the  stations  not  contained  in  the  third 
subset.  This  particular  way  of  gri  uping  stations  has  the  characteristic  that  each  subset 
samples  approximately  the  same  regions  of  the  Earth.  Therefore  it  is  expected  that 
the  bias  introduced  by  the  different  combination  of  stations  on  the  measured  eieen- 
periods  may  be  minimized.  The  fluctuations  in  the  eigenperiods  found  for  each  subset 
results  primarily  from  errors  due  to  noise,  data  handling,  and  an  unknown  elfcct  of 
combining  different  station  data.  The  excitation  criterion  was  applied  separately  to 
each  one  of  these  subsets  and  the  eigenperiods  were  measured.  In  general,  the 
difference  in  eigenperiods  among  subsets  is  very  small  and  indicates  that  the  excitation 
criterion  renders  stable  results.  In  many  cases  the  summation  of  only  20  spectra  gave 
stable  estimates  of  overtone  eigenperiods.  The  difference  between  the  average  eiecn- 
periods  over  the  four  subsets  and  the  eigenpenod  found  from  a  combination  of  all 
stations  is  given  in  Table  1.  The  standard  deviation  of  the  eigenperiods  determined 
using  all  stations  was  estimated  dividing  the  standard  deviation  found  for  the  subsets 
by  yJ2,  as  each  subset  contains  half  the  total  number  of  stations.  The  standard 
deviation  of  each  measurement  is  listed  in  Table  1.  Those  standard  deviations  are  of 
the  same  order  of  the  SEM  found  for  the  Alaskan  earthquake  free  oscillations 
(Dziewonski  &  Gilbert  1972a’)  for  those  modes  identified  in  both  cases.  This  agree¬ 
ment  indicates  that  summing  vectorially  the  spectra  or  taking  the  average  over 
measurements  at  single  stations  gives  results  of  a  similar  accuracy. 

The  average  standard  deviation  for0S, 3  to  0S5tl  eigenperiods  measured  at  different 
subsets  is  0-025  per  cent.  As  each  subset  samples  the  same  region  of  the  Earth,  the 
average  0-025  per  .ent  may  be  considered  the  minimum  expected  bias  introduced  by 
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determtniMoVeleCtl0n  °f  St3tlons  and  ePicentre  location  for  any  0S„  eigenpcriod 

Figs  29  and  30  show  the  0S„  and  0Tn  cigenperiods  measured  for  the  Alaskan  and 
reference*0  eVentS'  taking  the  avcraSe  cigenperiods  published  by  Derr  (1969)  as 

The  mean  difference  between  Derr's  cigenperiods  and  the  Colombian  data  in  the 
range  trom  0Sn  to  US4„  is  +0  030  per  cent  while  the  mean  difference  between  Derr's 
averages  and  the  Alaskan  data  is  0  058  per  cent.  The  Colombian  u  <ta  is  closer  to 
Derr  s  averages  than  the  Alaskan  data,  hut  while  the  lines  for  the  Alaskan  and  Colom¬ 
bian  shock  are  ronshlv  parallel  to  each  other,  they  show  a  more  erratic  departure  from 
Derr  s  averages  This  demonstrates  consistency  of  data  obtained  for  the  Alaskan 
and  Colombian  shock  and  also  suggests  that  the  Derr  s  data  have  lar^e  seattenne. 
The  systematic  dillercnce  between  the  eigenpcriods  on  the  Alaskan  and  Colombian 
shocks  in  the  range  from  nSn  to  is  0  045  percent.  It  is  smaller  than  the  standard 
deviation  ol  the  measurements  at  individual  stations  nut  it  is  larger  than  the  standard 

rTrt).0*  the  TCun  -ThlS  mdicatcs  tll:it  thc  heterogeneities  and  anisotropies  of  the 
fcarth  are  such  that  lor  a  given  epicentre  and  thc  present  distribution  of  stations  the 
average  eigentrequen.y  measured  over  all  stations  mav  differ  from  the  degenerate 
vc;“  °L fhe  " dcal  aycrage  laterally  homogeneous  Earth  bv  an  amount  larcer  than  the 

if  u  The,  U  circ  cs  m  Flt'  ”9  mdieate  tne  cigenperiods  measured  for  the  Colombian 
shock  on  the  co-latitudinal  component.  The  mean  deviation  between  the  vertical  and 
co-latitudmal  measurements  is  0  023  per  cent.  This  difference  can  be  attributed  to  the 
different  combination  ot  stations  used  in  each  case.  As  shown  in  Fie.  3  for  KOD 
the  excitation  ol  thc  vertical  and  horizontal  component  arc  different  tor  each  mode 
at  a  given  station.  The  horizontal  spectra  were  not  combined  with  the  vertical  spectra 
because  while  0S„  modes  dominate  the  vertical  component  they  are  more  contaminated 
by  other  modes  in  the  horizontal  component,  see  .-I,  and  |0  in  Fis>.  12.  If  we  had 
combined  both  data  thc  systematic  deviation  between  the  Colombian  and  Alaskan 

measurements  would  have  been  reduced  to  about  0  035  percent  but  it  still  would  be 
significant. 


-38- 


J.  A.  Mcndigurco 


ni '  "7  ~  °  f  0  ”  "!ea*urw  inc  Alaskan  and  Colombian  shocks 

potted  as  a  function  ol  w.  Tie  0  reference  line  corresponds  i0  ilie  average  eigen- 
periods  published  by  Deir  (I9'.9). 

11  is  remarkable  that  the  difference  in  eigen  periods  have  different  sien  for  0T„  and 

.°n  "ft,  1  Sr°  dlffcr^nt  m  m:|Pnitude.  The  systematic  difference  in  cTn  cieenperiods 

for  s  Z  frl°ni  "7?  'S  °:321  pcr  ccn{  or  7  ,in««  kir.ee r  than  the  difference 

or  0S„  cipcnptrioc.5.  flic  Colombian  data  shows  faster  Ray  lei  eh  waves  and  slower 

made  bv  A1  u  f  ,‘hC  ,Ai!,skan  s.hock’  A  sim,lar  »>'!«  of  observation  was 

waves  at  ?  u  „  ?  f°Und  a  d,fTeicncc  ,n  P:,3SC  'docity  of  0  4  per  cent  for  Love 

waves  at  -00  s  measured  over  many  great  circle  paths  for  the  Niieata  1964  event  and 
over  paths  throne, i  Pasadena  ( toksbz  &  Anderson  1963) 

anrf% 'V‘!h  °Ver,0RCS  l0r  A  modcs  arc  thcrc  suffeicnt  data  to  compare  Alaskan 
nbian  measurements.  The  Colombian  cieenperiods  ,S-  to  ,SM  are  con- 

0  non  ^  Sma  CP  i  tK  Alaskan  Thc  avcra£C  difference  is  0  115  percent  while 
0-090  per  cent  is  thc  average  S.D.  of  the  measurements. 


6.  Thc  excitation  criterion,  general  remarks 

caHv^h  TCeSS  °mthe  CXCltat,on  cntcrion  depends  on  how  well  we  can  predict  theoreti- 
2‘i  hv  0scllljtl0n  srpCCt,ra  at  3  Siven  statl011  Since  thc  theoretical  spectra  are 
derixed  b>  assuming  an  Earth  and  source  model  the  question  is  how  critical  is  the 
choice  ol  those  models  for  a  successful  peak  identification 

from" HB 1  ?'C  ' r,Mi'lla!TrS  ‘dcnt,licd  m  th,s  PaPcr  ‘he  change  of  model 

from  H B  6  (Haddon  ik  Bullcn  1969)  to  the  5  OS  A/  (Press  1968;  Kanamori  1970) 

m.lr°dUCC  any  s,?niticant  difference  in  thc  results.  The  main  features  of  the 
corctical  spectra  at  each  station  remained  thc  same,  and  thc  onlv  effect  of  chancing 
the  Earth  model  was  to  shift  the  cigcnpcriods  and  shghtiv  chance  the  amplitudes  This 

sha'rvl  ,Sf^C  ,0  thr lacl  th3t  3  Cl,anpc  in  Ear,h  model  does  not  chance  drastically  the 
shape  of  the  eigenfunctions,  at  least  for  modes  of  low  radial  order,  bor  higher  over- 
tones  the  shape  of  thc  cigcnlunctions  becomes  more  structure  dependent  and  in  those 
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SoBH^rT  .ldcnIlfica,io"  should  bc  expected  when  using  a  more  accurate  Earth 
assumed  Earth  model  CXC,,a‘'°n  Cri,Cri°n  is  not  c"‘*ca,,y  d^"d«’‘  ■"* 

A  more  serious  difficulty  is  posed  by  the  lateral  heteroeencities  of  the  real  Earth 
Mcndiguren  d972a)  observed  that  a  laterally  homogeneous  Earth  model  is  adequate 
t  predict  the  spectra  ol  inodes  with  lone  wavelengths.  But  for  shorter  wavelengths 
the  lateral  hetcro  enc.ties  cause  an  increasing  disagreement  between  theoretical  and 

r^ZZd  an  ,"'';!l,V  a‘  SOmc  Cr,tlcal  Mucicy  the  predicted  spectra  do  not 

resemble  the  observed  ones.  That  critical  wavelength  is  dnierent  for  different  modes, 
and  depends  on  the  uepth  ot  penetration  in  each  case.  E or  example  for  „  /'  modes  the 

K»T>Cax  Pien°d  'S  ai'°und  3ut)s  uh,lc  lor  o-S  modes  it  is  around  180  s  i  Mendiaurca 

72a).  In  general  those  modes  with  much  of  their  energy  in  the  upper  300  km  of  the 
Earth  arc  more  ntlcctcd  than  those  modes  penetrating  deeper.  I  her elore  for  each 
particular  mode  there  is  a  critical  frequency  beyond  which  it  is  not  possible  to  orcJict 
adequately  the  spectra  assuming  a  laterally  homogeneous  Earth  and  hence  the  excita¬ 
tion  criterion  cannot  he  applied. 

The  assumption  of  a  point  source  model  in  the  excitation  criterion  scheme  used 
or  the  Colombian  earthquake  was  adequate  because  the  linear  dimension  of  the 
source  were  smaller  than  100  km  i  Mcndiguren  1972a)  or  about  It)  times  smaller  than 
the  shortest  wavelength  ot  the  idcntiiied  mo  Ics.  For  shocks  associated  with  larger 
faul  dimensions  the  source  hmteness  should  be  taken  into  consideration  <e  e  Press 

of  a  .  I960)  when  computing  the  theoretical  spc.tra  and  the  phase  correction  before 
vectorial  summation. 

The  complications  introduced  by  the  presence  of  aftershocks  in  the  study  of  free 
oscillations  have  been  overestimated  in  the  past  (K.  Aki.  private  communication). 
Earlier  it  was  argued  that  as  the  amplitude  goes  like  the  square  root  of  the  energy  the 
amplitude  of  long  period  oscillations  should  be  proportional  to ./jo''-  Where  A/  is  the 
surface  wave  magnitude.  But  in  fact  the  amplitude  at  long  periods  is  proportional  to 
the  seismic  moment  (Aki  1967)  or  proportional  to  I o,/'. 

Anyhow,  when  applying  the  excitation  criterion  in  the  presence  of  la  rue  aftershocks 
their  effect  on  the  observed  spectra  could  be  taken  into  account  if  their  source 
mechanism  and  seismic  moment  were  known.  The  predicted  spectra  would  be  in  tins 
case  the  result  ol  summing  the  spectra  for  the  main  event  and  aftershocks.  Aftershocks 
were  not  observed  in  the  case  ol  the  Colombian  shock. 

An  obvious  extension  ot  the  excitation  criterion  would  be  to  sum  veetoriallv  the 
spectra  ol  many  earthquakes  recorded  at  many  stations  (Mcndiguren  1972c)  The 
spectra  for  each  shock  should  be  corrected  for  their  corresponding  source  mechanism, 
this  is  similar  to  the  super-resolution  technique  proposed  by  Gilbert  (1970)  to 
separate  singlets  of  each  multiple!. 

Records  of  different  components,  namely  horizontal  and  vertical,  or  records  from 
different  types  ot  instruments  could  also  be  combined.  It  should  be  possible  to  collect 
systematicaUy  data  lor  many  years  aid  continuously  improve  the  resolution  of  the 
method  and  the  accuracy  of  the  measured  eicenpenods.  flic  use  of  manv  earthquakes 
is  specially  desirable  because,  as  shown  in  Section  5.  the  epicentre  location  and  the 
arbitrary  distribute  ot  stations  introduce  a  bias  in  tiic  eitrenperiod  determination 
which  may  be  larger  than  the  standard  deviation  of  the  average. 

The  excitation  criterion  can  also  be  applied  to  other  types  ol  gcophvsical  observa- 
predicted^6  ^  CXCItatl°n  °‘  tlie  dltrercnt  superimposed  signals'can  be  theoretically 
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Appendix 

Data  processing 

Records  of  the  three  components  long  period  seismograph  at  82  WWSSN  stations 
were  digitized.  Since,  lor  most  ol  the  records,  the  shortest  period  present  was  longer 
than  8  s,  the  records  were  digitized  at  an  interval  of  4  s  to  avoid  the  aliasing  ctfect. 
In  a  few  cases,  where  high  frequencies  were  present,  the  digitizing  interval  was  reduced 
to  2s.  The  starting  time  of  digitization  was  taken  as  early  as  practically  possible,  but 
in  most  cases,  the  first  train  of  higher  modes  could  not  be  included  because  the  record 
trace  was  off  scale.  The  digitization  of  each  record  was  terminated  when  a  clear  train 
of  surface  waves  was  no  longer  visible.  The  resulting  record  length  ranged  from  4  to 
18 h.  The  digitized  records  were  corrected  for  recording  drum  drift  (Mitchell  & 
Landisman  1969).  In  order  to  save  computation  time  when  computing  the  spectra, 
the  number  of  digitized  points  was  reduced  to  one  fourth  by  successive  averaging 
over  9,  7,  5  and  3  points  (Blackman  &  Tukev  1958).  Then  one  out  of  every  four  points 
was  taken,  increasing  the  time  interval  between  points  to  16s.  In  addition,  the  digitized 
records  were  high-pass  filtered,  taking  out  those  periods  longer  than  701)  s. 

Since  the  theoretical  spectra  are  given  for  the  radial,  co-latitudinal  and  azimuthal 
components,  the  N-S  and  E-W  records  of  each  station  were  vectoriallv  summed  to 
synthesize  the  co-latitudinal  and  azimuthal  components.  The  re> ultant  records  were 
analysed  using  the  fast  Fourier  transform  over  16384  points.  In  the  next  step,  the 
spectra  were  corrected  for  instrument  response  (Hagiwara  1958)  and  for  the  low-pass 
filter  response  of  the  decimation  procedure. 

A  digitized  version  of  the  Nana  struinmetor  record,  with  readings  at  every  60s, 
was  provided  by  Drs  S.  Smith  and  D.  Anderson.  The  etTect  of  Earth  tides  was 
eliminated  through  high  pass  filtering. 
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All  the  theoretical  and  observed  spectra  shown  in  the  figures  of  this  study  were 
multiplied  by  ai*  for  an  adequate  graphical  representation. 

Comments  of  the  individual  peak  identifications 

The  measured  eigenperiods  were  considucd  reliable  and  reported  in  this  study 
only  when  the  measurement  was  not  allectcd  by  interfering  modes  or  noise.  Large 
peak  amplitude  compared  to  the  amplitude  of  the  surrounding  noise  and  a  smooth 
and  symmetrical  peak  shape  were  considered  indications  of  small  interference  In 
addition,  the  magnitude  of  the  standard  deviation  and  the  stability  of  the  measure¬ 
ments  among  dilfcrent  subsets  w  ere  taken  into  account  to  judge  the  reliability  of  each 
measured  cigenperiod.  An  additional  critical  test  belore  reporting  an  identification 
was  to  plot  the  measured  eigen  frequencies  as  a  function  of  n  (co-latitudinal  order 
number,  e  g.  0S„)  and  require  a  close  tit  to  a  smooth  curve.  All  these  elements  were 
taken  into  account  to  accept  or  reject  identifications. 

The  theoretical  cigenfrequency  for  the  HBI  6  Earth  model  corresponding  to  each 
mode  being  identified  is  indicated  in  Figs  14-28  with  a  larger  vertical  tick  on  the 
frequency  axis. 

0S  modes.  The  fundamental  spheroidal  mode  0S  dominates  the  vertical  displace¬ 
ment  spectra  above  160-s  period  and  the  identified  peaks  have  a  large  amplitude  over 
the  surrounding  noise. 

,S  modes.  These  modes  are  heavily  interfered  with  the  0S  modes.  The  reported 
peaks  are  those  which  satisfy  most  ol  the  tests  for  a  reliable  measurement.  For 
example  in  the  case  ot  tS}v,  Fig.  1 7,  the  identified  peak  is  surrounded  by  other  large 
amplitude  peaks  but  the  identification  is  reliable  because  the  measured  cigenfrequency 

is  exactly  the  same  for  each  of  the  four  stations  subsets,  see  Table  1. 

2S  modes.  In  general  these  are  modes  well  observed  on  the  co-latitudinal  com¬ 
ponent. 

3S  modes.  Spectral  peaks  have  been  observed  from  }S5  at  415  s  to  3Ssg  at  87  s. 
3S9  is  a  special  case  with  a  large  signal  to  noise  ratio  (Fig.  33)  but  a  large  standard 
deviation.  The  reported  cigenperiod  satisfies  the  required" smoothness  in" the  period 
i’.  n  plot.  Observed  on  the  co-latitudinal  component. 

*s  modes.  Observed  within  a  short  frequency  range,  its  peaks  are  well  developed 
on  the  co-latitudinal  component. 

jSo  modes.  These  modes  arc  well  excited  and  the  identifications  are  very  reliable. 
Observed  on  the  co-latitudinal  component. 

6S0  modes.  These  modes,  like  the  S  modes  were  observed  in  the  vertical  com¬ 
ponent,  and  they  show  more  interfering  noise  than  those  spheroidal  modes  observed 
in  the  co-latitudinal  component. 

0 T  modes.  Reliable  measurements  have  been  done  up  10  0T3().  For  larger  n  values 
the  lateral  heterogeneities  ot  the  Earth  allectcd  the  spectra  to  the  extent  that  the 
laterally  homogeneous  Earth  model  is  not  valid  any  more  (Mendiguren  1972a)  and 
the  excitation  criterion  fails  to  render  accurate  results. 

,T  modes.  Modes  up  to  ,  T,4  have  been  identified.  These  modes  contain  energy 
at  larger  depth  than  the  n  7*  modes  and  arc  not  so  aficcted  by  the  lateral  hetero°eneities 
of  the  upper  200-  300  km  (Meno.guren  1972a). 

2T  and  3  7  modes.  Only  a  few  isolated  groups  of  modes  have  been  observed  with 
sufficient  clarity  to  be  reported. 
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